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ABSTRACT:  This  is  a progress  report  on  an  investigation 

to  test  the  theory  of  Lewis  and  von  Elbe  (i)  for  the 
ignition  of  propellants. 

Information  from  two  types  of  experiment  was  utilized 
to  calculate  ignition  energies  per  unit  area  and  "linear" 
surface  temperatures.  At  the  University  of  Michigan  (2,3) 

M2  propellant  was  ignited  by  hot  gas  under  forced  convection, 
while  at  the  Franklin  Institute  (5,5)  ignition  occurred  by 
free  convection  (including  radiative)  heat  transfer.  From 
these  experiments  it  is  shown  that,  for  constant  gas  flow 
rate  (forced  convection)  and  constant  gas  pressure  (free 
and  radiative  convection),  the  ignition  energy  per  unit 
area  decreases  with  increasing  gas  temperature,  which  is 
in  accordance  with  the  theory.  The  ignition  times  are 
too  long  to  allow  determination  of  whether  a minimum 
ignition  energy  per  unit  area  exists.  It  is  proposed  to 
extend  the  investigation  to  the  short  ignition  times  pro- 
vided by  the  adiabatic  compression  technique  developed 
at  NOL. 
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the  amount  of  energy  per  unit  area  necessary  to  ignite  a 
propellant  under  given  conditions. 
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I.  TABLE  OF  SYMBOLS 


arbitrary  constant 

constant  pressure  heat  capacity 

diameter  of  propellant  grain 

energy  flux 

gas  mass  flow  rate 

reduced  heat  transfer  coefficient 

heat  transfer  coefficient 

thermal  conductivity  of  guo 

thermal  conductivity  of  propellant 

arbitrary  exponent- 

ignition  energy  per  unit  area 

time 

gas  temperature 
ambient  temperature 
surface  temperature 
assumed  surface  temperature 
temperature  above  ambient  &t  x and  t 
surface  temperature  above  ambient  at  t 
gas  temperature  above  ambient 
flctltuous  gas  "temperature” 
distance  from  surface  Into  propellant 
dlffustvlty  of  propellant 


dimensionless 
any  units 
any  units 
any  units 
Iba/hr  ft2 
1/om 

any  units 
any  units 
any  units 
dimensionless 
cal/cm2 
any  units 

n._ 

OK 

°K 

C‘K 

°K 

°K 

°K 

oR 

any  units 
any  units 
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(contd. ) 

r 

Ignition  time 

any  units 

/ f 

gas  viscosity 

any  units 

Subscripts: 

G "greatest" 
L "least" 
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II . INTRODUCTION 

The  ignition  theory  advanced  by  Lewis  and  von  Elbe 
states  that  the  amount  of  energy  per  unit  area,  q,  neces- 
sary to  ignite  a propellant  surface  depends  upon  the 
energy  source  temperature.  The  lower  that  temperature, 
the  greater  must  q be.  Also  as  the  temperature  Increases, 
q approaches  a minimum  value.  This  minimum  q is  thought 
to  be  approximately  equal  to  the  excess  enthalpy  associ- 
ated with  the  combustion  wave. 

The  purpose  of  this  investigation  has  been  to  test 
the  theory  and  to  demonstrate  the  conditions  under  which 
it  is  valid. 

The  equations  for  the  local  heat  transfer  coefficients 
for  forced  convection,  derived  by  the  University  of  Michi- 
gan (6),  enabled  us  to  test  this  theory  under  conditions 
which  are  similar  to  those  found  in  practical  ignition 
systems „ Because  of  the  type  of  dependence  of  these 
coefficients  upon  the  assumed  surface  temperature  of  the 
grain,  it  is  possible  to  calculate  "greatest''  and  "least" 
coefficients  for  every  experimental  run.  Hence,  for  every 
run,  this  permits  the  calculations  of  "greatest"  and  "least" 
value  of  q and  surface  temperature,  TB. 


III.  MATHEMATICAL  MODEL 

To  calculate  T8  and  q at  the  time  of  ignition, T , 
it  la  necessary  to  solve  the  classical  heat  conduction 
equation  for  a cylinder  being  heated  by  convection.  To 
simplify  the  problem  a little,  one  can  substitute  a semi- 
infinite  solid  for  a semi-infinite  cylinder.  This  is 
possible  because  the  diffusivity  of  the  class  of  solid 
propellants  la  extremely  small  and  the  order  of  ignition 
times  involved  does  not  allow  for  any  appreciable  tempera- 
ture rise  much  below  the  surface . Hence  * to  a temperature 
wave  propagating  into  the  solid  from  the  surface,  the  solid 
appears  vc  be  or  infinite  depth  i.e.  a » emi- infinite  solid. 
The  problem  is  reduced  to  that  of  finding  va , the  surface 
temperature  above  ambient,  and  a as  functions  of  t for  a 
semi-infinite  solid  being  heated  by  convection  and/or 
radiation  from  a medium  at  temp  V above  ambient;  the 
initial  temperature  of  the  solid  is  zero. 
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•«?T' 


The  equation  to  be  solved  Is 

(1) 


Bat  'd>*r 
b*  ‘ *Xx 


for  the  boundary  conditions 

+ A at  = V o±  x=o 

■av 

(2) 


3X 

at(x,0)  - 0 


The  solution  (8)  Is 


(3)  -V"  = erTC  ~ " * 


iM 


'rfc(jJS**VEi) 


where 


(°°  , * 

erf c | = -4-  I €~ 


Putting  x - 0 and  t «•  v , one  gets  for  the  surface 
temperature  rise  at  the  Ignition  time. 


(4)  AT(r)  = V ji-fc1*  trer4c  -Wicr  j 


To  calculate  q we  Integrate  the  flux  over  the  Ignition 
interval,  and  this  is 


(5)  ^ v>/*1  = j f C®,t)  cU 


where  =-K 


V»x 


We  have 


[O)  *x 


a-AVe*~t,M«rfc(fe*WS*) 


The  surface  flux  Is 

(7)  = KlV€irt  eric 
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Substituting  (7)  Into  (5)  yields 

= (Tfe  M 1 f,  it  (1VT- U A. 

. . > . > c YTT  l ! 

K -K.V  A. 


(8) 


J. 


J 


After  evaluating  the  second  Integral  by  parts  we 
finally  get 


(9) 


^ rT  «rfc  -I  | 


where  q will  replace  q(o,r  ) and  vs  will  replace  vs(  r) 
from  now  on. 


(9)  can  be  rewritten,  by  substituting  (4)  Into  It, 
giving  (9')  „ „ voM  £.  XvT?  - ^5.1 

o Xk  lVrr  V J 


IV.  EXPERIMENTAL  AND  CALCULATED  RESULTS 

It  Is  desired  to  evaluate  v3  and  q for  the  Michigan  and 
the  Franklin  Institute  data. 


The  Michigan  equation  for  the  heat  transfer  coefficient, 
H,  between  a moving  gas  stream  and  a right  circular  cylinder 
where  axis  is  normal  to  the  gas  flow  is 


uo> 


The  h used  is  eq(2)  is  defined  by 


(ii) 


A-  77 

rs 


Hence 


1 - 

— 


a x (c£\»  j&2\"  mV1 

k d w u / Vt;7 
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For  any  one  experimental  run  all  the  constants  of 
(12)  are  known  except  T's  which  increases  from  ambient 
tempera ture , T0.  Therefore  h depends  upon  our  knowledge 
of  T'a . But  by  eq(4),  we  must  know  h in  order  to  cal- 
culate Ts  - T0  + vB i This  does  not  make  the  situation 
hopeless,  for  by  the  following  method  we  can  demonstrate 
that  the  indeterminancy  in  v8  is  very  small  compared  to 
the  actual  vs. 

For  a given  r during  any  one  experiment,  the  values 
of  v8  and  q will  be  greatest  when  h is  greatest.  Similarly, 
the  values  of  va  and  q will  be  smallest  when  h is  least. 

From  (12)  it  is  easily  seen  that 

us) 

an^ 

0.!£ 

(no  al~  (4j) 

<* 

By  replacing  T'8l  by  T0  in  (13)  we  get  hQ  and  hence 
v8q  and  qG  by  (4)  and  (9).  Then  replacing  T'Bq  by  T0  + vSo 
in  (14)  we  obtain  hj,  and  hence  vg^  and  qj,.  We  now  have 

/vr%i  < <vrv  < at4<. 

'1Z>  ) %L  * % 4 


Analysis  of  the  data  that 


(16) 


« 'Vs 

t 


In  their  report,  the  Franklin  Institute  had  calculated 
the  surface  temperatures  for  the  Michigan  work  using  the 
latter' b equation  for  the  local  heat  transfer  coefficient. 
Since  their  values  of  v*  must  obey  (15)  and  (16),  we  have 
used  these  value#  in  («r)  to  calculate  the  q!s  which  also 
obey  (15)  and  (16).  For  the  range  investigated  qn  and  q^ 
differ  by  0.1  - 0.2  cal/cm2.  The  lower  limit  is  for  shorter 
ignition  times. 
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At  this  point  it  must  be  pointed  out  that  the  ignition 
time,  in  all  cases,  was  determined  by  actuation  of  a photocell 
by  the  propellant  flame. 

In  Table  1 are  given  the  calculated  q's  for  forced  con- 
vection using  an  80#  No  - 20 % Og  gas. 

In  Tables  2,  3 * and  4 are  given  the  calculated  q's  for 
free  convection  for  No*  and  air,  respectively.  The  v8's 
shown  are  those  previously  calculated  (5). 

In  Table  3 are  given  qo,  q^,  v8r,  vSL  and  the  experi- 
mental conditions  for  a different  set  of  forced  convection 
experiments  (3).  These  are  characterized  by  much  shorter 
ignition  times.  The  values  of  A - 0.8l  and  n » 0.54  were 
used  in  (10). 

Figures  1,  2,  3*  and  4 are  the  plots  of  data  in  Tables 
1*  2,  3,  and  4,  respectively. 

Tables  2,  3*  and  4 are  plotted  in  such  a manner  in 
Figures  5*  6,  and  7 that  each  Figure  represents  a constant 
pressure  level. 

The  physical  constants  used  are  reported  in  the  Michigan 
and  Franklin  Institute  reports. 


V.  INTERPRETATION  AND  CONCLUSIONS 

Since  the  model  used  allowed  for  energy  transfer  only 
from  a hot  gas  to  a propellant,  the  values  of  va  and  q 
represent  those  quantities  at  the  end  of  the  ignition  interval 
t , if  no  chemical  reaction  had  taken  place  in  the  propellant 

From  Figures  1,  2,  3*  and  4 and  Table  5*  one  can  see  that 
for  all  runs  at  constant  flow  rate  or  constant  pressure,  the 
values  of  q decrease  with  increasing  gas  temperature.  Here 
we  see  that  the  q necessary  to  ignite  the  propellant  surface 
must  depend  not  only  upon  the  energy  source  temperature,  as 
stated  in  Reference  1,  but  also  upon  the  flow  rats  or  pressure 
level  for  convective  processes. 

The  effect  of  pressure  upon  q varies  markedly  between 
the  three  gases  for  ?*aturai  convection.  In  Figures  2,  3* 
and  4 it  is  seen  that  if  we  consider  the  gas  at  constant 
temperature  an  increase  In  pressure  gives  rise  to  an  in- 
crease In  q for  helium  but  a decrease  in  q for  nitrogen  and 
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the  8o£  N2  - 20#  O2  mixture.  This  Inversion  cannot  be 
explained  at  the  present  moment.  We  can,  however,  specu- 
late as  to  the  cause. 

Suppose  propellant  ignition  was  determined,  in  addition 
to  the  assumed  convective  heat  transfer,  by  a gas-solid 
phase  chemical  reaction.  Since  ignition  times  are  deter- 
mined by  photocell  activation  due  to  the  propellant  flame 
radiation,  any  phenomenon  which  would  cause  the  delay  of 
the  flame  reaction  would  increase  r and,  by  (9)*  q.  The 
degree  of  diffusion  of  the  inert  gases,  He  and  No*  into 
the  gas  reaction  zone  next  to  the  propellant  surface  would 
vary  according  to  the  pressure  levels,  and  the  mobility  of 
the  inert  molecules  compared  to  those  of  the  gas  products 
at  any  particular  level.  Thus,  one  may  suppose  that  the 
diffusion  process  could  cause  the  inversion.  It  must  be 
admitted  that  this  is  a.  very  sketchy  description  of  the 
hypothesis . 

On  many  of  the  plots  in  Figures  2,  3*  and  4 the  300 
psig  curves  tend  to  fall  out  of  place  when  compared  with 
the  other  pressure  levels.  The  experiment era  (5)  feel  that 
due  to  experimental  difficulties  it  is  not  possible  to 
determine  with  certainty  whether  or  not  any  apparent 
exceptions  in  their  calculations  (hence  in  our's)  are 
significant. 

It  is  important  to  note  that  within  the  range  of 
experimental  results  used  the  lowest  value  of  q is  1.7 
cal/cm'2 , ‘this  value  is  in  agreement  with  minimum  ignition 
energies  obtained  from  previous  (9)  steady  state  burning 
results.  These  gave  a value  of  2.55  cal  cm2  for  the  excess 
enthalpy  associated  with  a combustion  wave  as  measured  by 
integrating  under  the  temperature-distance  curve  of  steady 
state  burning  of  nitrocellulose  strands.  The  value  of 
2.0  cal/cm2  for  the  energy  associated  with  the  combustion 
wave  was  obtained  ( 9 ) by  measuring  the  heat  drain  from 
a burning  powder  strand  to  a copper  strip  cemented  to  the 
bottom  of  the  strand.  However,  there  is  no  evidence  to 
indicate  that  if  the  gas  temperature  is  Increased,  the 
value  of  q at  ignition  could  not  conr.inue  to  decrease  with 
shorter  ignition  times. 


VI.  FUTURE  WORK 

At  NOL  (10)  we  have  been  igniting  various  propellants 
by  adiabatic  compression  since  this  method  bears  close 
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resemblence  to  practical  gun  and  rocket  ignition  conditions. 
We  can  obtain  ignition  of  equivalent  semi-infinite  pro- 
pellant cylinders  in  the  order  of  two  or  three  milliseconds. 
Attempts  are  being  made  to  calculate  the  ignition  energies 
and  a "linear"  surface  temperatures  for  such  very  short 
ignition  times  under  conditions  of  continuously  varying  gas 
temperature  and  pressure,  and  mass  flow  rate.  It  is  hoped 
that  these  conditions  will  allow  us  to  ascertain  whether 
there  exists  a critical  minimum  ignition  energy  per  unit 
area.  The  computations  are  necessarily  very  complex. 
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TABLE  1 

ignition  energies  per  unit  area  for 

M2  PROPELLANT  USING  605&N2-2OS&O2 

OBTAINED  FROM  UNIVERSITY  OF  MICHIOAN 
DATA  FOR  FORCED  CONVECTION 


5-0  SCFM 


V 

T 

h* 

V * 

vs 

q 

290 

7.65 

15.36 

185 

8.88 

391 

3.05 

16.79 

210 

6.23 

516 

1.14 

18.85 

19.04 

224 

4.01 

540 

0.768 

206 

3.00 

5&5 

0.767 

0.6c8 

19.87 

234 

3.25 

CryQ 

20.47 

234 

0 • v^c. 

646 

0.£35 

20.68 

232 

2.81 

712 

0.413 

21.45 

238 

2.55 

3,5  SCFM 

296 

9.25 

12,89 

1 fio 

m C.  0 

403 

3.10 

14.04 

197 

5.85 

529 

1 .26 

15.75 

211 

3.93 

635 

0.687 

17.08 

218 

2.95 

731 

0.447 

18.43 

233 

2.56 

764 

0.456 

18.42 

238 

2!  65 

2.0  SCFM 

322 

8.82 

9 .66 

ile 

8.63 

405 

4.47 

10.72 

6,73 

501 

2.07 

11.53 

191 

4.59 

632 

0.873 

12.70 

191 

2.93 

•Calculated  by  Franklin  Institute  (5) 
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TABLE  2 

IGNITION  ENERGIES  PER  UNIT  AREA  FOR  M2  PROPELLANT  USING  He 


OBTAINED  FROM  FRANKLIN  INSTITUTE  DATA  FOR  FREE  CONVECTION  AND 

RADIATION 


p 

V 

v** 

T 

h* 

v8* 

q 

0 

420 

599 

9.12 

3-57 

168 

8.38 

495 

767 

5.82 

3.61 

182 

7.21 

37° 

957 

3.52 

3.69 

188 

5.77 

645 

1186 

2.41 

3.72 

199 

5.12 

795 

1746 

1.08 

3.85 

213 

3.55 

50 

420 

532 

8.54 

5.73 

205 

9.88 

495 

665 

4.41 

5.78 

204 

7.04 

570 

815 

3.22 

5.82 

223 

£ cr» 

^ X 

645 

986 

2.02 

5.90 

227 

l3.28 

720 

1178 

1.28 

6.00 

228 

4.18 

795 

1406 

0.90 

6.00 

234 

3.61 

100 

395 

473 

8.20 

7.08 

206 

9.80 

4yp 

633 

4.30 

7.12 

223 

7 .72 

570 

768 

2.46 

7.20 

223 

5.79 

645 

922 

f 1.96 

7,24 

1247 

15.60 

11.20 

I203 

13,62 

720 

1095 

1.23 

7.32 

24p 

1*  l« 
H.Hf 

1 95 

0.92 

7.32 

258 

3.99 

200 

420 

490 

7.10 

9.16 

238 

10.67 

495 

601 

4.33 

9.24 

253 

8.75 

570 

724 

2.70 

9.26 

263 

6.95 

645 

861 

* i 

9.31 

266 

5.71 

720 

1014 

1.20 

9.34 

273 

h P '7 

~r  • | 

795 

1183 

0.86 

9.42 

282 

4.21 

300 

420 

479 

5.32 

10.83 

235 

9.20 

495 

585 

3.50 

10.90 

256 

7.92 

570 

700 

2.06 

11.00 

259 

6.14 

645 

828 

1.38 

10.96 

265 

5.13 

720 

969 

1.00 

11.02 

276 

4.55 

795 

1128 

0.67 

11.00 

275 

3.66 

♦Calculated  by  Franklin  Institute  (5) 

**V  Is  a "temperature”  which  includes  the  radiation  boundary 
condition.  To  calculate  vs  and  q.  V is  substituted  for  V 
in  equations  4 and  9.  See  Reference  (4)  for  details  on 
calculating  V, 
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TABLE  3 

IGNITION  ENERGIES  PER  UNIT  AREA 
FOR  M2  PROPELLANT  USING  N2 

OBTAINED  FROM  FRANKLIN  INSTITUTE 
DATA  FOR  FREE  CONVECTION  AND  RADIATION 


p 

V 

r 

h* 

vB* 

q 

0 

420 

895 

30.50 

1.35 

187 

16.77 

495 

1207 

15'.  40 

1.38 

192 

12.44 

570 

1588 

8.90 

1 .40 

202 

9.33 

645 

2062 

5.50 

1.42 

213 

7.68 

720 

2625 

3.70 

1.44 

228 

7.40 

795 

3306 

2.60 

1.46 

247 

6.40 

30 

420 

722 

17.50 

2.53 

200 

13.60 

495 

941 

9.30 

2.56 

205 

8.43 

570 

1201 

5.60 

2.60 

li 

8.12 

645 

3.70 

2.63 

227 

7.17 

720 

1886 

2,45 

2.65 

238 

5.86 

795 

2314 

1.80 

2.69 

257 

5,57 

10) 

420 

677 

12.30 

3.33 

202 

11.53 

495 

870 

6,80 

*5  Q7 
1 

208 

8.88 

->  * w 

i r\r\' 7 
i 

ll  1 A 
*T  # J.KJ 

It  «• 

215 

7.17 

645 

1369 

2.70 

3.42 

226 

5.91 

720 

1688 

1.85 

3.44 

237 

5.14 

795 

2053 

1 .40 

3.49 

257 

5.05 

200 

420 

631 

8.80 

4.52 

208 

10.17 

495 

799 

5.00 

4.59 

216 

7.89 

570 

995 

3.10 

4.62 

224 

6.36 

645 

1223 

2.10 

4.66 

235 

5.59 

720 

1488 

1.50 

4.68 

249 

4.99 

795 

1791 

1.10 

4.72 

264 

4.53 

300 

420 

615 

5.40 

5.43 

194 

7.44 

495 

771 

3.10 

5.55 

200 

5.81 

570 

952 

2.00 

5.60 

210 

4.74 

645 

1162 

1.35 

5.64 

218 

4,li 

720 

l4o6 

1.00 

5.64 

232 

3.77 

795 

1676 

0.75 

5.72 

247 

3.46 

♦Calculated  by  Franklin  Institute  (5) 
**See  note  in  Table  2 
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TABLE  4 

IGNITION  ENERGIES  PER  UNIT  AREA 
FOR  M2  PROPELLANT  USING  AIR 


OBTAINED  FROM  FRANKLIN  INSTITUTE 
DATA  FOR  FREE  CONVECTION  ANT)  RADIATION 


p 

V 

r 

h* 

v8* 

Q 

0 

420 

878 

25  = 80 

1 .40 

177 

14.38 

495 

1178 

11.84 

1.43 

173 

9.54 

570 

1547 

7.69 

1.46 

191 

8.34 

645 

2009 

5.22 

1.47 

211 

6.91 

720 

r'.r  1'  ' 

3.36 

1.50 

223 

5.39 

795 

3187 

2.17 

1.53 

231 

4.66 

100 

420 

664 

12.02 

3.50 

204 

11.57 

495 

851 

5.84 

3.56 

onn 

5-  W 

7.84 

570 

1073 

3.62 

3.57 

208 

6.39 

645 

1337 

2.41 

3.59 

219 

5.47 

720 

1646 

1.57 

3.60 

223 

4.59 

795 

2008 

1.25 

J.VC 

24y 

4.56 

300 

420 

607 

4.4o 

5.66 

184 

6.20 

495 

761 

2.94 

5-77 

200 

5.p4 

570 

9?9 

1.80 

5.81 

204 

£.40 

645 

1143 

1.36 

5,85 

00-3 

UU  4 

4.07 

720 

1377 

1.00 

5.89 

2?7 

3.  fo 

795 

1644 

0.74 

5.93 

249 

3.44 

^Calculated  by  Franklin  Institute  (5) 
”Ss£  note  lii  Table  2 
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